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Abstract 



When a high voltage is applied on an asymmetric capacitor, it experiences a force 
O acting toward its thinner electrode. This effect is called Brown-Biefeld effect 

(BB), after its discoverers Thomas- To wnsend Brown and Paul- Alfred Biefeld. 
Many theories have been proposed to explain this effect, and many speculations 
can be found on the net suggesting the effect is an antigravitation or a space 
warp effect. However, in the recent years, more an more researchers attribute 
' 5j ■ the BB effect to a unicharge ion wind. This work calculates the levitation force 

due to ion wind and presents experimental results which confirm the theoretical 
results. 
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The Brown-Biefeld (BB) effect has been discovered in 1920 by Thomas Townsend 
Brown and Paul Alfred Bicfcld during their experiments with Coolidge X- 
ray tube. They observed a thrust acting toward the thin electrode. Thomas 
Townsend Brown made an extensive research on this effect and wrote several 
patents 

There is a site dedicated to the BB effect Q and a site dedicated to experi- 
ments of this effect 5] , but very few theoretical works have been written_onJ;he 
subject. Some of those tried to explain the effect by electro-gravitation 
or thermodynamics 

However, in the recent years more and more sources 0> [ill H, Q attribute 
the effect to corona ionic air propulsion. In 13] a lot of experiments are de- 
scribed, and the thrust force is described by approximate formulas based on ion 
propulsion, and in llj a full calculation of the thrust, based on the jet of the 
corona ion wind is performed. The calculations in [ll| were very accurate, but 
seemed to require a substantial amount of computing power. 

In this work we also calculate the force due to ionic air pro puls ion, but we 
adopt a different approach, based on the Deutsch assumption [ljl, which has 
been extensively used in the unipolar charge flow literature 15 , 16 , [ijl [3 
[23, 24 j, but also criticized (see J. E. Jones et al. 20] and A. Bouziane et al. 



(21(), by showing big differences in magnitude and direction between the exact 
electric fields and the ones obtained with the Deutsch assumption. However our 
goal is to calculate the thrust force and for this purpose the Deutsch assumption 
proves useful. 

The Deutsch assumption states that the equipotential surfaces of the Lapla- 
cian problem are equipotential also for the Poissonian problem, only with dif- 
ferent values of potential. Of course under this assumption, the electric field 
lines of the Laplacian and Poissonian problems are in the same direction at any 
location. 

As we shall see, the calculations based on Deutsch assumption [ljj], within 
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the Warburg region |26ll27| result in formulas which fit very well the experiment. 

In Section 2 we explain the operation principle and the configuration which 
has been analyzed. We also present the equations that have to be solved. 

In Section 3 we present the basics of the Deutsch assumption. Sigmond 
made a profound analysis on the subject, and summarized the work of many 
researchers concerning the use of the Deutsch assumption. We will summarize 
those findings, and bring some highlights on this issue. 

In Section 4 we solve the Laplacian problem and calculate the capacitance. 
Those results are new, because as far as we know this Laplacian configuration 
has not been solved yet. Also, this Laplacian solution is needed for the Poisso- 
nian problem when Deutsch assumption is used. 

In Section 5 we solve the Poissonian problem using the Laplacian solution 
and Deutsch assumption and display the calculated force and current. 

In Section 6 we describe our experiment and show the obtained results. 

In Section 7 we compare the measured and calculated results and express 
our calculated results by approximate formulas. 

The work is ended with some concluding remarks. 

2. THE OPERATION PRINCIPLE AND THE CONFIG- 
URATION 



The picture of the lifter on which we did our experiments is shown in Fig- 
ure [TJ The lifter is based on a thin anode wire at positive high voltage, over a 
grounded flat cathode. As follows from the explanation below, for best propul- 
sion, the flat electrode must be vertical. We made some experiments with 
negative corona and the thrust is considerably lower for negative corona than 
for positive corona, fact which is confirmed by Blazelabs [13j. Hence we start 
the explanation for positive corona. 

The working principle is as follows: when a high enough positive voltage is 
supplied, pairs of positive ion and electron (which are always randomly created 
by incident photons) are accelerated: the positive ion toward the cathode and 



3 



Figure 1: (color online) The picture of the lifter we built and used for experiments. It consists 
of a thin anode wire over fiat vertical cathode and is of triangular form. 

the electron toward the anode. Some of those pairs recombine emitting a photon 
which by the photoelectric effect on a neutral atom, ionizes it, creating the 
electron avalanche. High energy electrons hitting neutral atoms, ionize them, 
contributing to the avalanche. This process arrives to a steady state consisting 
in a corona region around the anode (in which the ionization process happens) 
and in which there are positive ions and high speed electrons. Hence this region 
is neutral and typically narrow and outside it there is a unicharge positive 
ion drift moving toward the cathode. The positive ions transfer most of their 
momentum to neutral molecules and while the positive ions "feel" the force of 
the electric field and hence move toward the cathode (according to the positive 
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ion mobility coefficient) the neutral molecules keep the inertia of the momentum 
they received. If the flat cathode were horizontal (and infinite) , the jet of neutral 
molecules would hit the cathode and hence the forces on the anode and cathode 
would have been equal and opposite, hence no net thrust. 

However, the cathode being vertical, and given the fact that the ions transfer 
most of their momentum to neutral molecules, those do not hit the cathode, but 
form an air jet downstairs, and by momentum conservation the lifter senses a 
net force upwards. Hence, the thrust force is calculated as the total force on the 
space charge. 

Also we treat the whole region as unicharge, and neglect the corona thickness, 
which is small for the case of positive corona 
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About negative corona: negative corona starts like the positive corona by 
randomly created pairs of positive ions and electrons, only here the positive 
ions are accelerated toward the negative thin electrode (called here cathode) 
and the electrons are accelerated toward the flat electrode (called now anode). 
Also for negative corona some pairs recombine, emitting a photon, but this 
time the photoelectric effect created by it is on the thin wire surface, which 
being negative, easily releases electrons, creating the electron avalanche. Also 
electrons hit neutral atoms and ionize them, but because the electrons move 
outward the negative thin electrode, this happens at lower velocities than for 
positive corona, hence this part of the process is less dominant. Hence, the 
photoelectric effect on the thin wire surface being dominant, the negative corona 
is very sensitive to the ability of the thin wire to emit electrons. If the thin wire 
surface has irregularities only some parts of it emit electrons, creating tufts (see 
Peek [3], Fig. 77). Also for a limiting voltage, for which the thin wire surface 
is not able to steadily emit electrons, the process may be bursty. This happens 
because the positive ions attracted to the thin wire lower the strength of the 
electric field near it, and this field is restored only after the emitted electrons 
enter the flat positive conductor and arrive through the power source to the 
thin negative wire. If the process arrives to a steady state, the negative corona 
has two layers: the inner layer called ionization layer which consists of positive 
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ions and out flowing electrons emitted by the thin wire, the outer layer in which 
electrons flow out bouncing into neutral atoms and combine in negative ions. 
Outside those layers there is a unicharge flow of negative ions and electrons (see 
25| . Fig. 1). We suppose the main reasons for the levitating force being smaller 
with negative corona are the tufts (reflecting a nonuniform ability of the thin 
wire to emit electrons) and burstiness which may persist also after corona builds 
up. We will not further deal with negative corona and all the results obtained 
in this paper are valid for positive corona only. 

Now the shape of the lifter must not necessarily be triangular, it may be 
rectangular or any other shape. Neglecting edge effects, the operation is de- 
scribed by a thin anode wire, over a vertical conducting plane. Hence we deal 
with a two dimensional problem described in Figure [2j 



T 2a/ 
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Corona wire 



Flat cathode electrode 



Figure 2: The simplified 2D configuration. The corona wire is at coordinates (0, 0) and its 
radius is a. The distance between the electrodes is called b — a, and the cathode length is c. 



Our lifter has the following dimensions: a — 0.075mm, b = 2.8cm and 
c = 4cm (see Figure [2]). Each side of the triangle is 0.2m, hence if we calculate 
the thrust force per unit length, we have to multiply by the perimeter of 3 x 0.2 = 
0.6m to find the total force. 

As long as the potential difference is below the corona inception voltage [3| , 
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there is no space charge, hence we have a Laplacian problem, defined by 



V 2 V L = 0, (1) 

where the index "L" denotes the Laplacian solution. The boundary condi- 
tions are Vl — Vq on the anode wire surface and Vj, = on the cathode wire 
surface, where Vq is the applied voltage. The electric field is El — — VVj,. 

In presence of space charge, different kinds of ions have different mobilities, 
defining the velocity of each ion as its mobility times the electric field. However, 
we use the average mobility [3, [ijl Q] for positive air ions, known to be p 
2 x 10~ 4 m 2 jv sec. Diffusion can usually be neglected 
non diffusive drift of ions is described by 



3, 



29, 30], hence the unipolar 



J = pv = npE Pl (2) 

where p is the space charge per unit of volume and Ep is not the Laplacian 
field, but the Poissonian field, influenced by the space charge itself via Gauss 
law: 



e V -E p = p. (3) 

Being a stationary problem, the current conservation condition V • J = 
must hold. So the Poissonian problem may be formulated by: V • {pEp) = or 



V • (V 2 Vp Wp) = 0. (4) 

Clearly, this equation contains the 3rd derivative, hence one needs an addi- 
tional boundary condition, and what is usually assumed is Kaptzov hypothesis 

ail- 

Physically Kaptzov hypothesis means that once the potential difference has 
been raised sufficiently for the corona to start, the electric field near the corona 
conductor remains constant and equals the inception value, even when further 
raising the potential difference. 
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Our goal is to calculate the thrust force, i.e. the x component of the total 
force on the space charge JJJ Ep x pd T, where Ep x is the x component of Ep, 
and the integral is on the whole free space. Using eq. ©, the total thrust force 
may be expressed as 



/ - / / / •/,<! '. {-,) 



and for this goal we need a full solution of eq. (U). Felici 



29( analyzed 



equation ((4]) and gave some particular solutions, and Feng 32[ calculated an 
exact solution of eq. ((4]) for concentric cylinders. 

As mentioned before, eq. (j4} is easier to solve by using the Deutsch assump- 
tion Q and this item is discussed in the following section. 

3. ANALYSIS OF THE DEUTSCH ASSUMPTION 



The Deutsch [14J assumption (DA) states that the equipotential surfaces of 
the Laplacian problem, defined by eq. ([T]) are equipotential also for the Poisso- 
nian problem, defined by eq. (j4|), only with different values of potential. This 
assumption is not valid in general, hence its usage violates some laws of physics. 
Sigmond Q studied the consequences of using DA, and showed that there 
are two general approaches in using DA. One approach called T-type, keeps 
V x Ep = 0, but violates current conservation. The other approach called L- 
type, satisfies the current conservation V • J — but violates the orthogonality 
between flux lines Ep and potential surfaces (i.e. results in V x Ep ^ 0). 

For more details, the reader is referred to [ijj], formulas (13)-(39), but for 
convenience we show here the main highlights of the issue. Dealing with a two 
dimensional problem we shall use orthogonal coordinates u, v describing the 
field direction and equipotential surfaces, respectively (the equivalents of I, a in 
[ill). Accordingly, the unit vectors will be noted as u, v. 

The assumption that the equipotential surfaces of the Laplacian problem 
remain equipotential in the Poissonian problem implies that the Poissonian field 
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Ep can be expressed as 



E P = 8E L , (6) 

where 6 is a scalar function. It is crucial for our solution to have a correct 
current distribution, because from eq. ([5]) we see that only a correct current can 
result in a correct force, hence we use the L-type described above. Requiring 
current conservation and using eq. @ and ^ we obtain 



V • J = V • (w9E L ) = V(>p0) ■ E L + (npO)V ■ E L = V(jip0) ■ E L = 0, (7) 

because V • El = 0. Now El pointing in the u direction, and being non 
zero, implies: 

^p6) = 0, (8) 
which means that [ipd must be a function of v only: 

(j,p6 = K{v). (9) 
Using again eq. ^ and © we get 

J = K(v)E L . (10) 



22] and he showed that one 



This method has been carried out by Popkov 
may obtain the potential difference on each field line as a function of the current 
density on the collector, and used cross-plot technique to obtain the current 
density in terms of the applied voltage. Sarma and Janischewskyj 2s| followed 
the same technique as Popkov, but used an iterative procedure to calculate the 
current density on each field line in terms of the applied voltage. A. Bouziane et 
al. 



24] ) used the above methods to evaluate the current density and field profile 
for several geometries, and found the results in agreement with experiment. 
Our goal is to calculate the thrust force in eq. (JSJ) , hence we will not follow 
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the above techniques, but rather look for an adequate function K(v). It will 
be helpful for example if we could choose this function to be constant in some 
regions of v. 

Now because we used the L-type, V x Ep ^ in general, so if we require 
V x Ep = we end with a contradiction. 

As we shall see in the next section, the Laplacian field lines are almost 
straight in the paraxial region around y — (see Figure^. It is therefore of 
interest to examine what happens in a region where the Laplacian field lines El 
are straight. In such region the unit vector u does not change with u, so is a 
function v only. Knowing that V x El = 0, results in El being a function of u 
only. 

We shall see that in this case we may require V x Ep = 0, because this 
belongs to the cases described in [lj3]> eqn. (38) and (39), for which the DA 
holds exactly. So requiring: 



= Vx£ P = Vx (9E L ) = V0 x E L + 9 V x E L = V0 x E L (11) 



results in V# being in the u direction, hence 9 being only a function of u. 
Therefore Ep — 9El is also a function of u only, and the same is true for 
p = e V • E P . 

The outcome is that pp9 in eq. © which is only a function of v must also 
be only a function of u. This is possible only if the function K(v) in eq. (|10p is 
a constant. 

This is helpful, since we may start our Poissonian calculations (Sect. 5) in 
the paraxial region with K(v) constant, and let it drop to when approaching 



the Warburg 



26, 



27] limit region. 



This is in accordance with the experimental knowledge that the current drops 



to outside the Warburg region 



34 



35 



36 



37| and also with Ieta 16 . who 



showed that the Laplacian solution reconstructs well the Warburg distribution 
for a pin to plane geometry. As we shall see, this method results in a calculated 
thrust force which fits experiment. 



10 



So, for using Deutsch assumption in the Poissonian problem, we need the 
solution to the Laplacian problem first, and this is done in the next section. 

4. SOLUTION OF THE LAPLACIAN PROBLEM 

There are many ways to solve the Laplacian problem (analytic or numeric) 
and we chose to solve it analytically, by separation of variables in cylindrical 
coordinates. 

The schematic configuration in Figure [2] does not allow to properly define 
boundary conditions in separate variables in cylindrical coordinates, hence we 
will make a slight change to this configuration. 

Anyhow the configuration in Figure [2] is not accurate, because it does not 
show the width of the flat cathode. In practice the cathode must have a finite 
width, and more than that, the curvature radius of the cathode at the location 
(x = b,y = 0) ( see Figure [2]) must be big enough so that no corona can be 
formed there [28| . 

The modified configuration is described in Figure 




X 

Figure 3: Configuration adapted to cylindrical coordinates. The cathode contour is described 
by the lines r = b, and <p = ±a. 
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The flat cathode surface is now defined along ip — ±a, where a is a fixed 
small angle of about 1.3° (based on the thickness of the cathode in our lifter). 
This model still includes edges at (r = b,p> — ±a), but those edges do not exist 
in the real engine and except of them the model is quite close to reality. At 
last, we take the cathode length c to infinity: the main interaction is between 
the electrodes and using a finite length c would still be solvable analytically, but 
would be an unnecessary complication. 

We define the potential V\ for the region a < r < b and < ip < 2n and the 
potential V 2 for the region r > b and a < p < 2tt — a. The boundary conditions 
are: 

V 1 (r = a,<p) = V , (12) 

where Vo is the applied voltage. Because of the mirror symmetry around the 
x axis, we require 

V 1 {r,<p) = V 1 {r,2TT-<p). (13) 
The potential continuity at r = b gives 



Vi(r = b,ip) = < 



< (p < a 

2tt - a < p < 2tt (14) 

V2(r = b, ip) a < ip < 2n — a 
and the normal field continuity at r = b requires 

d r Vi(r = b,a < tp < 2n- a) = d r V 2 (r = b, a < p < 2tt - a). (15) 
The potential is on the sides of the big conductor 

V 2 (r,(p = a) = V 2 (r,tp = 2ir-a) =0 (16) 
and must go to at r — > 00 : 
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V 2 (r oo,a < ip < 2?r - a) = 0. (17) 

We will use the well known solutions for the Laplace equation in cylindrical 
coordinates, for the z independent case, given by the trivial solution D — Ehir 
(where D and E are constants) plus the non trivial solution: 

^2{A v r" + B v r- v )(C v cosM + D v sin(^)), (18) 

V 

where one may consider only non negative values of v : because negative v 
just switches the roles of A v and B„. 

Let us start with V 2 - Condition (fT7| excludes the trivial solution and the 
r v solution which diverge at r — ¥ 00. Also, the requirement V2(r, ip = a) = in 
cq. (fl6l) imposes a combination between the sin and cos terms in eq. (fT8|) of the 
form sin(i^(( ( 9 — a)) and the constant B v may be normalized for convenience to 
B„/b~ u . Hence we may write the following expression for V%: 

V 2 = Y, B Ar/b)-"sm(is{<p-a)). (19) 

To satisfy the requirement V^ir, if = 2t: — a) = in eq. (|16[) . we need 
sin(i/(27r — a — a)) = 0, or v(2tt — 2a) = run (where m is a positive integer), 
thus giving the values of v = y i_a/ 7r • So the expression for V2 may be written 
as: 

li-fw)^*,^). (20) 

Now we look for an expression for Vi . Because the non trivial solution is ip 
dependent for any r, to satisfy condition (|12j) . we need also the trivial solution. 
Also, because V\ is defined in a region with circular continuity, adding 2ti to 
(p must result in the same potential, hence v must be an integer, say m. For 
satisfying eq. (|13l) only the cos solution must be taken, and again, we may 
normalize the constants so that the power is on r/a instead of r, so we may 
write the following expression for V\ : 
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V^D-Elnr+Y^ ( A m{r/a) m + C m {r/a)~ m ) cos(to^). (21) 

To satisfy condition (|12j) , the to ^ terms of the non trivial part of eq. ([2Tj) , 
must be identically for r = a. Given that the functions cos(to^) are orthogonal 
in the interval < <p < 2tt, each term of the series must vanish for r — a, hence 
we get A m = —C m , for to 7^ 0. The m = term gives just a constant, which 
may be absorbed in the trivial solution, but it proves convenient to name the 
to = term Lq, and to scale separately the trivial solution so that it results 
in an arbitrary constant V for r = a, and for r = b, so that D — Ehir = 

V (l - §03) ■ So we obtain for Vi: 

Vi =V'(l- ^tM) +Lo+jTAm ((r/o)"* - (r/a)— ) cos(to^). (22) 

\ \ I ) / m—l 

It is to be mentioned that we do not lose any generality with the above 
scaling of the trivial solution, because after scaling, the trivial solution plus Lq 
result in V + Lq for r — a and Lq for r = b, and the relation between them has 
not been established yet. 

The reason for choosing this approach is that all the unknowns, namely A m 
and Lq in eq. (|22j) and B m in eq. (|20[) must be proportional to the applied 
voltage Vo, so we may calculate them with the aid of an arbitrary V, and 
obtain V\(r = a,(p) = V + Lq, which can be scaled eventually by a factor 
Vo/(V + Lq) to be equal to Vo- So we may set from now ox\ V' = 1, and 
remember to multiply everything by Vq/{\ + Lq). One may verify that the 
alternative approach of absorbing Lq in the trivial solution results in much 
more complicated equations. 

Here one has to require the boundary conditions and solve 2 sets of matrix 
equations to find the vectors A mi B m and Lq. The details are worked out in 
the Appendix, and Figure |4] shows the potential in a 3D plot, for the physical 
values of our lifter, where Vq is normalized to 1. Because of the mirror symmetry 
around the x axis, we drew the potential only for positive y. 
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Cathode 



x [cm] 



Figure 4: 3D plot for the calculated potential. The applied voltage Vo is normalized to 1. 
The calculations used the physical values of our lifter, i.e. corona wire radius a = 0.075mm, 
distance between electrodes b = 2.8cm and the angle a = 1.3°. 



We may derive now the electric field. For region 1 we obtain the radial held: 



E 9Vl = V ° 1 

Or 1 + Lor 

and the circular held: 



1 



ln(6/a) 



mA m ((r/a) m + (r/a)-" 1 ) cos(m^) 

(23) 



m— 1 



Ei v = --^ = £ mA m ((r/a) m (r/a)- 1 ) sin(m^). (24) 

r oip 1 + Lq r 
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For region 2 we obtain the radial field: 



dV 2 Vp 1^ 1-1/2 ■ ( n , m ^- a 

E 7r = — — = - — -> —Bi(r/b) sin (i - 1/2) 

r t — ' 1 — ct/n V 



<9r 1 + Lo r 1 — a/n \ 1 — a/nJ 

(25) 

and the circular field: 



«,..!«. «if izi/^, (r/6) -S£ cos („ _ 1/2) . 



<p — a 



a/n 
(26) 

It would be useful to calculate the capacitance. The charge per unit of 
surface on the anode is r/ = E$Ei r {r — a, if), and the charge per unit of length 
is given by A = a d<p rj. Clearly, the integral on <p zeroes the sum in eq. (|23|) 
and we are left with: 

27re a V a 

m(6/a)(l + £ )' 1 ' 

So the capacitance per unit of length is 

r' = — = 27re ° (2*\ 
v Hb/ a )(i + L y [a> 

Of course, for a = it, from eq. (|A.8[) we have Lq = and we recover the 
known formula for the capacitance per unit of length for concentric cylinders. 

We are of course interested in small a, so we calculated the values of Lo for 
different ratios b/a, and fit an approximate formula for it - see Figure [5j 

It comes out that for a < 7r/100, Lq can be expressed as 

1.3035 , N 

L ^lnW^ +0 - 011 ' ^ 
resulting in the following approximate formula for the capacitance per unit 

of length: 



~ 1.011 ln(6/a) + 1.3035' v ; 
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Figure 5: (color online) Values of Lq for a = 1.3°. The stars are the exact calculated values, 
and the continuous line represents the fitted curve values given by eq. 1291 This fitted curve is 
valid for a < tt/100 rd = 1.8°. 



If we set the values for our lifter: a = 0.075mm and b = 2.8cm into eq. (|30|) 
we get C — 7.63pF/m, and multiplying by the lifter's perimeter 0.6m, the 
calculated capacitance comes out Calculated — 4.57pF. To check the validity 
of this result we measured the capacitance of our lifter and got C mea surod = 
4pF ± 5%, which is quite close to the calculated result. 

It is also useful to get a relation between the applied voltage Vq and the field 
intensity on the anode. Of course, the field on the anode wire is not constant, 
and depends on (p. But for a thin anode (b 3> a), the field is almost constant 
(see also discussion in the next section). 

One may verify that for r — a the absolute value of the sum in eq. (|23|) is 
much smaller than ln ^ a ^ hence we may write the approximate expression: 
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Va 



(1 + L Q ) aln(b/a) « a(1.3035 + 1.011 ln(6/a)), 



(31) 



E Xr (r = a, ip) 

where the second expression is a further approximation which uses eq. f)29[) . 
Again, for a = n, L$ = 0, and we recover the known expression of this 
relation for concentric cylinders. 



In Figure |6] we show the equipotential surfaces and the field lines - those are 
needed in the next section. 






x [cm] 



4 f 6 

Cathode 



Figure 6: (color online) The thin lines emerging from the anode are the electric field lines 
while the thick lines are the equipotential surfaces. The coordinates u and v denote the local 
direction of the field and the equipotential lines, respectively. The dashed line represents 
the electric field line which delimits the Warburg region, hence passes through the point 
(x = b,y = b/ tan(60°)). Only the y > region is shown, because of the mirror symmetry 
around the x axis. 



The Warburg 



2G 



region can be seen also in FigureEJ and will be referred 



to in the next section. 



5. THE POISSONIAN PROBLEM 
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In this section we use the Laplacian results obtained in the previous section 
to solve the Poissonian problem, i.e. the state of the system when the applied 
voltage is bigger than the corona inception voltage. 

F.W. Peek [28[ made an extensive research on corona inception for different 
geometries like concentric cylinders, parallel wires, etc. and published the results 
in his book. For all the configurations involving corona around a thin wire of 
radius a, the electric field on the surface of the wire at which corona begins (at 
room temperature) is given by B, = 3 x 10 6 (1 +p/y/a) V/m, where p ss 0.03 y/m 
with very slight variations of about 2% for different geometries, with different 



asymmetries for the electric field (as found by Peek [28J, page 63). 

We do not have the exact value for our geometry, but as explained by Peek 
himself, before corona starts, the field very close to a thin wire behaves like the 
field on the surface (which is almost constant if the wire is thin) times the wire 
radius, divided by the distance from the wire. So when having the above Ei 
value on the wire surface, one can easily find that at a distance of py/a from 
the wire surface the field is 3 x 10 6 V/m, this way assuring a field of more than 
3 x 10 6 V/m in a wide enough region around the wire to allow corona to start. 

So we can safely use Peek formula: 

£, = 3 x 10 6 (1 + 0.03/v^) V/m (32) 

as has been done a lot in the literature, for different configurations of thin 
wire electrode near any other electrode [ljlH, 32 1. 



Given a — 0.075mm, we know that for our case Ei = 13.392 MV/m. 

Now we can find at which voltage the corona starts, i.e. the Corona Inception 
Voltage (CIV). The CIV is the voltage for which the Laplacian field on the 
surface of the corona wire equals to Ei. For this we do not have to use Peek 
formula for CIV, we have the Laplacian solution for our problem. Using the 
approximation (|31j) . results in 7.32 KV or running the solution and measuring 
the exact relation results in 7.42 KV, so the difference is less than 2%. Hence 
we may use: 
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CIV = 7.4 KV. 



(33) 



Because of the asymmetry around the anode, the field on the corona surface 
is not completely uniform, and that is the reason for the slight variations in 
p in Peek formula for different asymmetric configurations. But fortunately, 
the corona wire being very thin, the field on the corona wire surface is almost 
uniform (up to variations of 0.5%), so we do not have to worry about it. (See 
also discussion before eq. (|3"Tj) ). 

So for a voltage bigger than the CIV, the Laplacian solution is not valid 
anymore, and we need the solution to the Poissonian problem. This requires 
the simultaneous solution of eq. ([3]) and ©, where for J we use eq. (fTU|) . 

The coefficient K in eq. (|10l) is unknown, but given the fact that the Deutsch 
assumption is accurate for y = 0, i.e. where the field lines are straight (sec 
discussion in Sect. 3 and eqn. (38) and (39) in Sigmond [18j). one can iterate 
the coefficient K in this region to get J Edx = Vq, where Vb is the potential 
difference for which we solve. 

The numerical solution is described in Figure [3 which is a zoom on a region 
of Figure [5] 



Figure 7: (color online) Zoom on a region Figure \E\ The calculations are done along the 
Laplacian field lines, and start with n = on the anode surface up the final value of n on the 
cathode surface. The lengths on the sides of each area are called dO, dl, L0 and LI. 
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The "P" prefix for Poissonian is omitted for brevity, and having a two di- 
mensional problem J represents the current per unit of perpendicular length, 
and p represents the charge per unit of surface. The iterative calculation is done 
between each pair of field lines (see Figure [7]) by: 

Pn = J n /(jiE n ), (34) 
where J n is J at the location of E n , and it is known for a given K, and 



E n +l — {PndavLav/tO + dO E n )/dl, 



(35) 



where d av = (dO + dl)/2 and L av = (LO + Ll)/2. Those are the numerical 
implementations of eq. (J2J) and ©• The initial condition for the iteration is 



Kaptzov assumption 



31] (explained in the introduction), which requires: 



E =E l . (36) 

After finishing the calculation between the first two field lines (i.e. in the 
region of very small y we check the result of J^E dx E n L aVn . Say its value 

n 

is 2Vq, we have to reduce K by a factor of approximately 4 (approximate because 
the initial condition for E is independent on K). This process converges very 
quickly (3-4 iterations), and after establishing K we can process the calculations. 
As explained in Sect. 3, we let K drop to when reaching the Warburg 
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limit region (see Figure [6]) . Within each area element we also calculate 



the x component of the force 



F Xn = E XnPn , (37) 

where E Xn is the x component of E n . The total force is eventually summed 
on the whole area, and knowing J we sum across the line fields (in the v direc- 
tion), obtaining the total current I. 

In the final stage the force and the current are multiplied by 2 to account 
for the symmetric y < region and the values being per unit of length of the 
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lifter, are multiplied by the perimeter 0.6m. The force is normalized to show 
the lifted mass in grams. 

The values of Vq for which we did the calculations, have been chosen to 
correspond to the values on which the experiment has been done (see next 
section). The calculated results are presented in Table [T] 



Table 1: Calculated results 



Vq [KV] 


Current [mA] 


Mass that can be lifted [g] 


11.12 


0.084 


1.54 


12.9 


0.145 


2.65 


14.08 


0.19 


3.55 


14.4 


0.219 


4 


15.64 


0.27 


4.98 


16.8 


0.34 


6.16 


18.9 


0.46 


8.42 


20.5 


0.586 


10.74 


21.8 


0.678 


12.44 


22.8 


0.76 


13.9 



6. THE EXPERIMENT 



The diagram of the experimental setup is shown in Figure O 



High 
voltage 
source 




R 




Lifter 



Figure 8: The diagram of the experimental setup. The lifter is connected to a high voltage 
source through a resistor. For different values of resistors, the input and lifter voltages are 
measured. The lifting force is measured in the cases it was big enough to lift. 

We used the lifter shown in Figure [1] for our experiments. The corona wire 
is a regular copper wire of radius 0.075mm, and the cathode is an aluminum foil 
of width 4cm (called a and c in Figure [21 respectively). The sticks that hold the 
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device are of balsa wood, and the perimeter of the engine is 60cm. The above 
materials are relatively light, so that the mass of the lifter is 7g. 

For different values of resistors, the input and lifter voltages are measured 
and the lifting force is measured. The lifting force has been measured only for 
the cases for which the lifter lifted, by counter balancing it, knowing that its 
mass is 7g. 

The results of the measurements are shown Tabled] The current is calculated 
using the 2 measured voltages and the resistor. 



Table 2: Measured results (NA means not available) 



R [MQ] 


Lifter voltage V Q [KV] 


Source voltage [KV] 


Current [mA] 


Lifted mass [g] 


308 


11.12 


25.6 


0.047 


NA 


154 


12.9 


25.58 


0.082 


NA 


110 


14.08 


25.51 


0.104 


NA 


88 


14.4 


25.46 


0.126 


NA 


66 


15.64 


25.42 


0.148 


NA 


44 


16.8 


25.3 


0.193 


NA 


22 


18.9 


25 


0.277 


7.5 


10 


20.5 


25 


0.450 


10 


6.8 


21.8 


24.92 


0.459 


12.2 


3.3 


22.8 


24.8 


0.606 


13 



7. COMPARISON AND APPROXIMATED FORMULAS 



One may see that the calculated forces fit well to the measured forces, with 
an average deviation of about 6.5%. This suggests that our method gives a good 
estimate for the lifting force. 

Let us first analyze the relation between force and current, by comparing it 
with a simpler case of straight, parallel and uniform field lines. If the field lines 
are in the x direction, our coordinates u, V (see Figure [6]) correspond to x, y 
respectively. For this case, applying eq. <j5j> F — JJJ J x d 3 r and integrating 
over the y and z coordinates yields F = — J dx. The integral just gives the 
length of the field lines, let us call it d, hence the force on the ionic space charge 
is F = (Id/fi). 

Certainly, if this configuration is a parallel plates capacitor, the above force 
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on the space charge does not produce thrust, because the wind hits one plate. 
But one may build a lifter with approximate straight, parallel and uniform 
field lines by using many thin wires on top of many flat vertical aluminum foil 
cathodes (see 13(). And for such a lifter F = (Id/ fi) describes well the thrust 
force. 

Examining the connection between our calculated force and current in Ta- 
ble [T] we find that the relation between force and current is 18.4 g/mA ± 0.8%, 
or in terms of force instead of mass, we get 179.72 N/A. As explained above 
for parallel field, we would expect this relation to be proportional to 6, which is 
the distance between the electrodes in our lifter. 

Hence we repeated the procedure explained in Sect. 5 for different values of 
6, and summarized the results (which include the original value of 6) in Table [3l 



Table 3: Relation F/I for different distances between electrodes 



Distance relative to original 


Calculated F/I 


F/I approximated by 1.2846/ fj, 


0.5 


89.42 


89.88 


0.75 


134.41 


134.82 


1 


179.72 


179.76 


1.5 


269.39 


269.64 


2 


358.82 


359.5 



Obviously, one can express this relation with an excellent accuracy by: 

7(1.2846) 



F 



/' 



(38) 



Now we analyze the current- voltage relation for the calculated and measured 
results, which are shown in Figure [SI 

Both measured and calculated results may be curve fitted with the formula 



I = KiVoiVo - CIV), 



(39) 



although the measured results look farther from their fitted curve than the 
calculated results (might be because of some technical problems during experi- 
ment). 
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Figure 9: (color online) The current-voltage curve of the lifter. The '*' describe the measured 
results and the '+' describe the calculated results. The curves are the closest approximations 
to the above results using / = K\ Vo(Vb — CIV). For the calculated results K\ = 2.16 mA/KV 2 
and CIV = 7.42KV, while for the measured results K 1 = 1.82 mA/KV 2 and CIV = 9.46KV. 



The curve fitted to the calculated results used K\ — 2.16 mA/KV and 
CIV = 7.42KV, while this fitted to the measured results used K x = 1.82 mA/KV 2 
and CIV = 9.46KV. 

We remark that the measured results show a too high CIV and this is con- 
nected to the fact that the correspondence between measured and calculated 
results is poor for low currents. This might be explained by the fact that the 
CIV is affected by additional factors like temperature, air pressure, or partial 
damage to corona wire, because of repeated operation |28| . On the other hand, 
higher values of current and voltage are less affected by the CIV (see eq. (f39| ) , 
and there we get a much better correspondence between measured and calcu- 
lated current. 

For a point to horizontal plane Warburg [26] found the constant K\ to be 
proportional to the distance between the electrodes to the power —3.17. This 
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power has been corrected to —3 by Jones 38}. 

For concentric cylinders, Townsend [33J found the constant K\ to be pro- 
portional to the capacitance, and the distance between electrodes at power —2, 
namely: K\ — 4/iC/R 2 , where fi is the ion mobility and R is the outer cylinder 
radius. 

We would expect our case, being two dimensional, to behave similar to 
Townsend's concentric cylinders, i.e. K\ should be some constant multiplied 
by ixC/b 2 . This assumption is tested with the results obtained for different 
values of b, after curve fitting, and summarized in Table [4] 



Table 4: Constant Ki of the relation / = KxVoiVo - CIV) (in units of mA/KV 2 ) for different 
distances between electrodes 



Distance relative to original 


Curve fitted K\ 


K\ approximated by 1.866/1(7/6^ 


0.5 


9.07 


9.63 


0.75 


3.89 


4.04 


1 


2.16 


2.17 


1.5 


0.96 


0.91 


2 


0.53 


0.5 



We see that K\ can be expressed by the approximate formula 



K x = 1.866fiC/b 2 (40) 
with an accuracy of about 5%. 

We may combine the eqs. ([38 ]) . (|39 ]l . ([iO]). with eq. (|3l) and (|32j) to get 
an approximate formula for the thrust on the lifter as function of the applied 
voltage Vq: 



F ^ 4J927re 4 ( i.oii m ( J!°) + 1.3035 - 3 x 106 (° + -° 3 ^) < (41) 

where all the magnitudes are in MKS units. 

We remark that the mobility /i canceled out. This does not mean that 
eq. (|4ip works for negative ions, because as explained in Sect. 2, the calcula- 
tions done in this paper are valid only for positive corona. Also we see that the 
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force increases as the distance b decreases. 

CONCLUSIONS 

In this paper we calculated the force on a levitation unit, built as a thin 
anode wire over a vertical cathode plane, as the electric force on the space 
charge and compared with experiment. 

Our calculations were based on the Laplacian solution for this electrostatic 
configuration, which by itself is a new result. 

With the aid of the Laplacian solution and Deutsch assumption we were able 
to calculate the force on the levitation unit and the current. The calculated 
results showed a good correspondence with the measured results. 

Also, we derived for this type of lifter the relation between thrust and current 
in eq. (|38p and the current-voltage characteristic curve in equations (|39p and 

601). 

And finally based on the above relations, we derived an approximate formula 
for the force as function of the applied voltage in eq. (|4Tj) . 
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Appendix A. The calculation of the Laplacian solution 

First we require condition (ITS]) : 



i 00 00 _ m R 



■ sm 



m <£> — a 
i/tt) \ 2 1 — 01/71", 
(A.l) 

and the above condition holds for a < ip < 2ix — a. In this range, the 
sin (yT^oJt^) functions are orthogonal, so we multiply eq. (|A.1|) by sin 
(for any positive integer n) and integrate on ip over the range [a, 2n — a]. We 
use the following integrals: 



n ip — a 







sin — )d(p = < (A.2) 

2 l-a/Kj I 1(1 - a /7r) n odd 



, \21 — a/7r/ \21 — a/iT J 

where 5 mn is the Kronecker delta. We also use: 



2-7T — Ct 

n ip — a 



sin n 1 r~ cos(mip)dip = G nm , (A.4) 

1 2 1 — a/TT / 



where G nm is a matrix defined by 



n even 

cos(ma) (( » /2) "( 1 1 _~%^) ) 2_ m 2 m^(n/2)/(l-a/n) n odd . 
— (7T — a) sin(ma) m = (n/2)/(l — a/ir) n odd 

(A.5) 

The last case defined by m = (n/2)(l — a/n) is not likely to happen if a — >■ 0, 
except for very specific values of a, but we calculated this case for completeness. 
It is to be mentioned that for this case a/n = 1 — n/(2m), for some specific m 
and n, hence ma — mn — nn/2 is equivalent to the points tt/2 or 3ir/2 on the 
unity circle, so that cos ma = 0, and this has been used in the above calculation. 
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Also, sin(ma) could be written as (— i)™-™/ 2 - 1 / 2 . 

After performing the above integrals, we obtain from eq. (|A.1[) the following 
result for the B coefficients in eq. (|20p : 



1 ] n even 

-nnBr, 



2 ' 4 



(! - a A)RW + ^-=1 G «™ mA ™ ((V«) m + (V«)^ m ) n odd 

(A.6) 



This result implies that B n are for even n, and this is expected because of 
the mirror symmetry around the x axis. 
Now we require conditions (114[) : 



< if < a 

L a+T,~ =1 A m {(b/a) m -(b/a)- m ) cos(rn V )= ^ Q 2tT - a < if < 2vT 

£m=l g»nsin(^ fgi^r ) Q! < <y9 < 27T - a 

(A.7) 

The above condition is defined for < ip < 2tt. In this range, the cos(my) 
functions are orthogonal, so we multiply eq. (|A.7[) by cos(ny) and integrate on 
ip over the range [0, 2tt], for n being any non negative integer. 

First, if we do it for n = 0, the sum on the left side vanishes, and using 
eq. (|A.2[) . with m instead of n, we obtain the following expression for Lq: 

L = -(l-a/n) ]T % (A.8) 

m—l (odd) 

and we could have even omitted the word "odd" , because we already know 
that the even B coefficients from eq. (|20l) are 0. 

Now we multiply eq. (IA.7[) by cos(nip) and integrate from ip = to 2ir, for 
n =/= 0. This time L$ vanishes, and after using cos(rw^) cos{rrup)dip = 7r5 mn , 
and the result (|A.4|) for m and n switched, we obtain from eq. (|A.7j) the following 
result for the A coefficients in eq. 
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irA n ((b/a) n - (b/a)- n ) = ]T G mn B m . (A.9) 

m=l (odd) 

We redefine now the odd indexes (like m in eq. IA.9[) by 21 — 1, I going from 1 
to oo. Hence the vector B m and the matrix G mn (like in eq. IA.9j) are redefined 
to Bi and G lm . 

This applies to equations ([20]). (|A.6j) . (IA.8I) and (|A.9|) . 

We rewrite now eq. (|A.6[) and (|A.9|) in matrix form, to solve 2 matrix 
equations with 2 unknown vectors for the coefficients A and B, then find Lq 
with eq. (|A.8|) and eventually calculate the potential and scale it by the factor 
V /(l + L ). 

Hence we define the diagonal matrix Q, by its components as Q m n = 
{{b/a) m - (b/a)~ m ) S mn , and obtain: 



ttQA = g T B (A.10) 

and we define the diagonal matrices V, by P mn = {{b/a) m + (b/a)~ m ) S mn , 
S, by S mn = (2m — l)5 mn and J , by J rnn = mS mn . We also define the column 
vector K by its components K m = 1, obtaining: 



- \kSB = -4(1 - a/7r)—^--S- 1 K + QJVA. (A.ll) 
2 ln(o/aj 



We isolate B from eq. (|A.11[) : 



B= -(l-a/vr)— \—S- 2 K--S- x QJVA (A.12) 
and set it in eq. (IA.10[) to obtain a closed form solution for A: 



a = -(i - 1— (ttQ + -g T s- 1 gjvr 1 g T s- 2 K. (A.13) 

7r ln(o/aJ 7r 

Clearly, for a — ir (concentric cylinders), A, B and Lq are all 0, hence V2 = 
and we are left only with the trivial solution for V\, as expected. 
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